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We have developed a near-field microwave microscope to locally excite a superconducting film 
and measure second and third harmonic responses at microwave frequencies. We study the local 
nonlinear response of a YBa2Cu307_f thin film grown on a bi-crystal SrTiOa substrate. The 
location of the bi-crystal grain boundary is clearly identified by the microscope through higher 
harmonic response, and the spatial resolution is on the order of the magnetic loop probe size, 
about 500/im. The harmonic power and spatial resolution are successfully modeled with a one- 
dimensional extended Josephson junction simulation. From the model, the 2nd order harmonic 
response is dominated by Josephson vortex generation and flow. A geometry-free nonlinear scaling 
current density Jnl — W 4 ~ 10 5 A /cm 2 is also extracted from the data, indicating that the grain 
boundary weak link is the dominant nonlinear source in this case. 

PACS numbers: PACS: 74.25. Fy, 74.25.Nf, 74.62.Dh, 74.72.Bk 



The nonlinear properties of high-T c superconductors 
have been of great concern in microwave applications, 
although the microscopic origins of the nonlinear re- 
sponse still remain uncertain. [Q All superconductors 
have an intrinsic nonlinearity associated with the non- 
linear Meissner effect (NLM E) . Calculations of harmonic 
response based on BCS theory,]^, || Ginzburg-Landau 
(GL) theory, Jl and microwave field-induced modulation 
of the super/normal fluid densities near T c 0] have been 
proposed to describe the intrinsic nonlinearities of su- 
perconductors. Extrinsic sources of nonlinearity include 
grain boundaries, || edge effects, and weakly-coupled 
grains. Q Many experiments have studied the intermod- 
ulation power, || [| harmonic generation, fiot |ll) or the 
nonlinear surface impedance of superconductors jl^, 
as a function of applied microwave power. However, 
most nonlinear experiments are done with resonant tech- 
niques, which by their nature study the averaged nonlin- 
ear response from the whole sample rather than locally. 
Such techniques usually have difficulty in either avoid- 
ing edge effects, which give undesired vortex entry, or 
in identifying the microscopic nonlinear sources. There- 
fore, a technique capable of locally measuring nonlinear 
properties of samples is necessary for understanding the 
physics of different nonlinear mechanisms. In addition, 
most existing experimental techniques focus on 3rd or- 
der nonlinearities, which can be conveniently studied by 
intermodulation techniques, but rarely address the 2nd 
order nonlinear response. Here we present a technique 
to locally characterize 2nd and 3rd order nonlinearities 
through spatially localized harmonic generation. 

In prior work,^] we studied the intermodulation sig- 
nal from a high-T c superconducting microwave resonator 
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using a scanned electric field pick-up probe. Both the 
"global" and the "local" intermodulation power mea- 
sured with the open-end coaxial probe were presented. 
However, the local measurements were actually a su- 
perposition of nonlinear responses that were generated 
locally but propagated throughout the microstrip and 
formed a resonant standing-wave pattern. To avoid 
this loss of spatial information, we have developed a 
non-resonant near-field microwave microscope, to non- 
destructively measure the local harmonic generation from 
un-patterned samples. This technique has the additional 
advantage of operating equally well through T c and into 
the normal state of the sample. 

In our experiment (Fig. [I]) , low pass filters are used to 
filter out higher harmonics generated by the microwave 
source, and guarantee that only the selected fundamental 
frequency is sent to the sample. A loop probe, which is 
made of a coaxial cable with its inner conductor forming 
a semi-circular loop shorted with the outer conductor, 
induces microwave frequency currents of a controlled ge- 
ometry in the sample. The probe can be translated over 
the surface of the sample in the x-y plane. The har- 
monic signals generated in the sample couples back to the 
loop probe. An advantage of this probe is the localized 
and directional microwave current induced on the sample 
surface, which enables the study of direction-dependent 
nonlinearities. Reflected harmonics are selected by 2 high 
pass filters, before being amplified by ~ 65 dB, and mea- 
sured by a spectrum analyzer. 

The sample is a 5001 thick YBaaCuaOy-,? (YBCO) 
thin film deposited by pulsed laser deposition on a bi- 
crystal SrTiOa substrate with a 30° mis-orientation. The 
distance between the loop probe and the sample is fixed 
by a I2.5/im thick Teflon sheet placed between them. 
Measurements of the temperature dependent 3rd order 
harmonic power (P3/)are first performed, both above 
the grain boundary (GB) and far away from the grain 
boundary (non-GB), as shown in Fig. 0. A strong peak 
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in P 3 f(T) is observed at T c = 88.9K (measured by ac 
susceptibility) at all locations on the sample. The peaks 
have similar magnitudes at both locations although there 
is a slight shift of T c . Note that all measurements are 
taken near the middle of the film where we have veri- 
fied that current-enhancement edge effects are absent. Jf4| 
This peak near T c is predicted by all models of intrinsic 
nonlinearities of superconductors, e.g. the NLME,[|[ |j| 
GL theory, Q etc., and a power-3 dependence of the P3/ 
on the input microwave power (-P/) is observed, in agree- 
ment with the predictions of those models. The order pa- 
rameter of the fragile superconducting state is easily sup- 
pressed by the microwave current at temperatures near 
T c , so that its nonlinear dependence on the external rf 
current dominates the microwave response. Also note 
that no feature is seen in P2/(T) near T c , as expected for 
a time-reversal symmetric superconductor. 

The SrTi03 substrate is a nonlinear dielectric at low 
temperatures, and we have measured harmonic re- 
sponse from bare SrTi03 substrates below 80K.|Q The 
nonlinear response is confined to narrow temperature 
ranges at temperatures when the substrate becomes res- 
onant due to its temperature-dependent, high dielec- 
tric constant. However, at temperatures below 80K, a 
strongly temperature dependent P 3 j is observed above 
the YBCO bi-crystal grain boundary, while no detectable 
P3/ is seen away from the grain boundary. Such clear 
contrast is evidence that the observed P3/ is from the 
grain boundary, not the SrTi03 substrate. Power de- 
pendencies of P2/ and P3/ were also performed at both 
locations at 60K (< T c ) and 95K (> T c ). As shown in 
Fig. H, strongly power dependent P2/ and P3/ are only 
observed above the grain boundary, while no response 
is seen above the background noise level away from the 
bi-crystal grain boundary. 

To demonstrate how the microwave microscope is able 
to spatially resolve a localized source of nonlinearity, a 
measurement of P2/ and P3/ along a line crossing the 
grain boundary is performed. As shown in Fig. ^, a clear 
peak in both P2/ and P3/ is observed above the GB, 
with a width of about 500/im . The width of the observed 
-P2//-P3/ peaks are about the size of the loop probe, which 
determines the spatial distribution of the surface current 
on the sample. This interpretation is confirmed by repro- 
ducing this peak with the extended resistively shunted 
Josephson junction model (ERSJ) discussed below. 

It is well known that applying microwave current to 
a single resistively shunted Josephson junction gener- 
ates harmonics at all odd integer multiples of the drive 
frequency. || ^| To obtain a more comprehensive under- 
standing of weak link junctions, the ERSJ model was 
introduced to model long Josephson junctions, such as 
the YBCO bi-crystal grain boundary. || 0, [l2| In this 
paper, we present ERSJ models to simulate a YBCO bi- 
crystal grain boundary as either an array of identical in- 
ductively coupled, or independent (uncoupled), Joseph- 
son junctions acting in parallel. From prior work with 
these junctions, we expect each to have a critical cur- 



rent of 8/j.A, and a shunt resistance of 20f2, and the 
currents applied to each junction vary according to the 
surface current distribution on the film induced by the 
loop probe. We sum up the nonlinear potential differ- 
ences across all junctions, and extract the higher har- 
monics from this collective nonlinear potential difference 
via Fourier expansion. || The spatial distribution of the 
surface current density is calculated from a simplified an- 
alytical model of an ideal circular loop in a vertical plane, 
with radius 27Q^Lm, coupling to a perfectly conducting 
horizontal plane 382.5/im away from the center of the 
loop. The magnitude of the current density is deter- 
mined by a much more sophisticated microwave simula- 
tion by the Ansoft High Frequency Structure Simulator 
software, which also produces a similar surface current 
distribution. 

The calculation of the uncoupled ERSJ model is per- 
formed by Mathematica, and is shown as the dashed line 
centered on 4mm in Fig. ||. This model predicts a narrow 
distribution of P3/ of greater magnitude (almost 20dB) 
than is observed experimentally. It also predicts no P2/, 
due to the absence of Josephson vortices in this model. A 
power dependence calculation from this uncoupled ERSJ 
model is also performed and compared with experimen- 
tal results (dashed line in Fig. ||). The comparison shows 
qualitative agreement with the P3/(P/) data. Calcula- 
tions with the inductively coupled ERSJ model, which 
includes Josephson vortices, performed by WRSpice, give 
very good descriptions in both magnitude and spatial res- 
olution of the experimental results for both P2 / and P3 f 
(solid lines in Fig. [I|). 

To further evaluate the capability of our microscope 
to detect intrinsic superconducting nonlinearities due to 
different mechanisms, we extract a geometry-free scaling 
current density, Jatl,|| from our data. We assume that 
the quadratic nonlinearity in the kinetic inductance of 
our sample dominates the nonlinear response through 

Lki = ^ [ X L (T,J) 2 J 2 dV = L Q + ALI 2 , (1) 

•*■ J volume 

and P 3/ = ( " A 4 £j3 ) 2 /(2Z ), where P 3/ is the 3rd har- 
monic power in the sample, u> and I are the microwave 
frequency and current, Zq is the characteristic impedance 
of the transmission line, Lki is the kinetic inductance 
of the sample, Lq and AL are the linear and 2nd or- 
der nonlinear terms in Lki, Al(T, J) = Al(T, 0)(1 + 
( JwUT) ) 2 ) 1 ^ 2 ^ s t ne temperature and current dependent 
London penetration depth, and Jnl is the nonlinear- 
ity scaling current density. Different microscopic mod- 
els of nonlinearity predict different values and tempera- 
ture dependences of Jnl{T)- For example, in the nonlin- 
ear Meissner effect and Ginzburg-Landau theory, Jnl ~ 
10 8 A I cm 2 or higher, ji], || except for temperatures close 
to T c , while the Jnl of a long 1-D Josephson junc- 
tion array in series combination is expected to be about 
10 6 Aj cm 2 or less.Q In our case, we follow Booth's 
algorithm Wq to calculate the nonlinear term in the in- 
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ductance, and the expected 3rd order harmonic response, 
all as a function of Jjvx> using the current distribution 
calculated from the model of the loop probe. However, 
it should be noted that the calculation of J ml is based 
on the assumption that the nonlinearity in the kinetic in- 
ductance of the superconductor is purely quadratic. As 
long as P3/ does not scale with Pj?, the calculated Jnl 
will be power dependent.^ 

Assuming X L {T = 0, J = 0) = 1500A, we obtain the 
Jnl calculated near the grain boundary to be Jnl ~ 
1.5 x 10 5 A/cm 2 , while the sensitivity of our setup is cur- 
rently limited to J NL < 1.3 x 10 6 A/cm 2 .[^4) However, 
it is suggested from the model calculation that thinner 
films and stronger coupling between the film and the loop 
probe will give stronger nonlinear response from the same 
mechanism, and improve the sensitivity to the nonlincar- 
ities associated with larger values of Jnl- 

We have created a novel near-field microwave micro- 
scope that is capable of locally measuring the nonlinear 
microwave behavior, thus providing an important new 
tool to study nonlinearities of superconductors. The mi- 



croscope is able to distinguish GB and non-GB areas 
through spatially resolved nonlinearity measurements of 
higher harmonic generation. The spatial resolution is de- 
termined by the size of the loop probe, which is about 
500/iTO in the current setup, and can be improved by 
reducing the size of the loop probe, and decreasing the 
loop/sample separation. A quantitative understanding of 
the magnitude and spatial resolution of the 2nd and 3rd 
order harmonic response is well achieved by the coupled 
ERSJ model. A geometry-free scaling current density 
Jnl is extracted from our measurement, and can be used 
to further study the physics of different nonlinearities of 
superconductors. 
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FIG. 1: Experimental schematic. The 6.5 GHz microwave 
signal is low-pass filtered before being sent to the sample, 
and high-pass filtered before being amplified and measured 
by the spectrum analyzer. 



FIG. 2: Temperature dependence of 3rd (a) and 2nd (b) 
harmonic generation on non-GB YBCO surface (circles) and 
YBCO grain boundary (GB) (triangles) in the temperature 
range 5~90 K. The fundamental frequency is 6.5 GHz and 
input power +8dBm. 



FIG. 3: Power dependence of 2nd and 3rd harmonic gen- 
eration on both non-GB YBCO surface and grain bound- 
ary at 60K (<C T c ) (filled-triangles: P 3 f on GB; diamonds: 
Pif on GB; circles: P3/ on non-GB; invered-traingles: P2/ 
on non-GB). The uncoupled-ERSJ model calculation for the 
power dependence of the 3rd harmonic generation on the grain 
boundary at 60K is also plotted as a dashed line for compar- 
ison. 



FIG. 4: Second- (-P2/, circles) and third-harmonic (P3/, tri- 
angles) power along a line crossing the bi-crystal grain bound- 
ary, which is located at X=4mm. Calculations of P 2 / and P 3 / 
from the coupled ERSJ model are shown in solid lines, and 
P.if from the uncoupled ERSJ in a dashed line. Inset is the 
schematic of the probe-sample arrangement. 
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Fig. 3 S.C. Lee, et al. 
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Fig. 4 S.C. Lee, et al. 



